Three unequal masses on a ring and soft triangular billiards 
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The dynamics of three soft interacting particles on a ring is shown to correspond to the motion of 
one particle inside a soft triangular billiard. The dynamics inside the soft billiard depends only on 
the masses ratio between particles and softness ratio of the particles interaction. The transition from 
soft to hard interaction can be appropriately explored using potentials for which the corresponding 
equations of motion are well defined in the hard wall limit. Numerical examples are shown for the 
soft Toda-like interaction and the error function. 
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Hard walls in billiard systems were extensively 
modelled by (^-potentials since they allow for sim- 
ple analytical relations of velocities and angles be- 
fore and after the collisions with the walls. How- 
ever, equations of motion are not well defined at 
the collision point. To analyze the transition to 
soft walls, which are more realistic, it is essen- 
tial to have well defined equations of motion since 
in general no simple analytical solutions are ob- 
tained. The present work suggests that appropri- 
ated soft walls are those for which the forces, not 
the potentials, become (5-functions in the limit of 
hard walls. This allows for better numerical in- 
vestigation of the soft-hard transition. A general 
scaled Hamiltonian is derived for three unequal 
masses interacting particles on a frictionless ring, 
which nicely describes this transition and shows 
that the dynamics occurs inside a soft triangle 
billiard. 



I. INTRODUCTION 

Although physically realizable potentials are inher- 
ently soft, theoretical models usually focus on hard po- 
tentials since they allow for analytical solutions. Billiard 
systems with hard walls are one example of such theo- 
retical models and have been extensively study and well 
understood. For example, it is well know [l]-[3| that the 
motion of three particles on a frictionless ring with point- 
like interactions is equivalent to the motion of one par- 
ticle moving freely inside the triangular billiard which 
hard walls. In such description the hard-walls potentials 
are represented by 5-functions which makes the collisions 
with the walls very simple, and analytical results can be 
obtained for the real dynamics and also in the tangent 
space. However, not much has been done in the descrip- 
tion of the transition from hard to soft interactions. Soft 
walls apparently do not destroy trajectories found in the 
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hard-wall limit Q and induce the appearance of regular 
islands in phase space HhS]. Some works about the dy- 
namics in soft billiard have to be mentioned, one particle 
motion in an atom-optic billiard [13] , quantum conduc- 
tance in the soft wall microwave billiard lljf . the effect 
of wall roughness in granular Poiseuille flow [12| and how 
the confinement of the equilibrium hard-sphere fluid to 
restrictive one- and two-dimensional channels with soft 
interacting walls modifies its structure, dynamics, and 
entropy [l3l [Tij and nonlinear dynamics in general 15 1 . 

There is a key issue for an appropriate description of 
soft to hard transitions in the context of billiards. Usu- 
ally the hard walls are modelled by <5-functions which 
allow for a simple description of the collision processes, 
since variables (angle and velocity) before and after the 
collisions with the walls can be given explicitly. However, 
the corresponding equations of motion are not well de- 
fined. In the case of soft walls the variables before and af- 
ter the collisions cannot be given in general explicitly, and 
equations of motion must be solved numerically. Thus, 
in order to study continuously the transition from soft to 
hard walls, equations of motions must be well defined in 
all cases. A first contribution Q in this direction ana- 
lyzed the dynamics of two interacting particles inside a 
ID billiard with soft walls, where the soft walls were mod- 
elled by the error function, which describes continuously 
and correctly the mentioned transition. The purpose of 
the present work is to generalize these results to the case 
of three soft interacting particles on a ring, and analyze 
the correct transition to point-like collisions. A general 
scaled Hamiltonian is derived which nicely describes this 
transition and also shows that the dynamics occurs in- 
side a soft triangle billiard. The dynamics depends only 
on the masses ratio between particles and softness ratio 
of the interaction. 

The paper is organized as follows. While in Section ITT1 
the hard wall triangle billiard is revisited, Section IIIII 
presents the general soft Hamiltonian system for three in- 
teracting particles, our main analytical result. In Section 
II VI we apply the soft Hamiltonian to the Toda model with 
unequal masses. Section [V] presents the example of two 
suitable potential which could be used to appropriately 
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describe the soft to hard transition for three particles. 
Finally Section [VII summarizes our main results. 



II. REVISITING THE HARD TRIANGULAR 
BILLIARD 

It is well know [J-[3| that the motion of three particles 
on a frictionlcss ring with point-like interactions is equi- 
valent to one particle moving freely inside the triangular 
billiard with angles 
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and colliding elastically with the sides of the triangle. 
Here M = m\ + 7712 + 7773 and m, is the mass of the 
particles (i = 1,2,3). The Hamiltonian can be written 
as 



H B = K M + 5(q 2 - qi) + S(q 3 - q 2 ) + S(qi -q 3 + L). (1) 

where Km = v\l^ m \ + P 2 /^ m 2 + is the mass 

dependent kinetic energy. The collisions occur at q± = 
<?2 j qi = <?3 and qi = q 3 + L where L is the circumference 
of the ring. The point-like collision between particles 1 
and 2 defines one side of the triangle at qi — q 2 = 0, 
and the collision of these particles with particle 3 defines 
the other two sides of the same triangle. For 7713 — > 00 
(?/ = 7r/2) we get the right triangular billiard which corre- 
sponds to the motion of two particles mi and 7772 moving 
inside the ID box with hard walls. In this case the inter- 
action between particles 1 and 2 is the point-like collision 
and the fixed particle 3 plays the role of a ID hard-wall. 
In such systems the Lyapunov exponent is zero (l6l . [l7j 
and the whole dynamics can be monitored by chan ging 
the angles of the triangle billiard pj. It was shown |l8j| 
that the Yukawa interaction between particles 1 and 2 
is enough to generate positive Lyapunov exponents. It 
is also worth to mention it is possible to relate the 
linear instability inside the triangular billiard with the 
Lyapunov exponents from quadratic irrational numbers, 
which arc related to the angles of the triangle, and thus 
to the masses ratio. 



III. THE SOFT TRIANGULAR BILLIARD 

As observed in Hamiltonian (TTJ), the collisions with the 
fixed hard-wall can be represented by 5-potentials. How- 
ever, the corresponding equations of motion are not well 
defined. Therefore, to describe analytically the transi- 
tion to hard walls, we include soft interactions between 
particles which, in a given limit, are expected to describe 
the collisions with the hard-walls. Let us start with the 
Hamiltonian of the three particles on a frictionlcss ring 
given by 



H B = K M + V 12 
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where er^ = 1, 2, 3) defines the softness of each pair- 
wise interaction. It is assumed that the interaction po- 
tential Vij between the particles depends only on the re- 



lative position between them and that for o~i 



a 



5-likc function is obtained for the corresponding force, 
and not for the potential. In addition, bounded motion 
is expected below a certain energy which will be specified 
later. Before going into details about the appropriated 
potential Vij which could be used, we rewrite Hamilto- 
nian @. Using the orthogonal transformation 
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for the three particles on a ring, and the linear transfor- 
mations x' = fix, y' = j3y, z' = /3z, dr/dt = /3 = , 
the final scaled Hamiltonian reads (without primes) 



K + V 12 (y) + V 23 (x,y) + V 31 (x,y), (4) 



where K = p 2 x /2 + p 2 y /2 + p 2 z /2 is the mass independent 
kinetic energy and 



V 12 (y) = V 12 (-y). 



V 23 (x,y) = V 23 
V 3 i(x,y) = V 31 
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Here /7 2 i = m 1 m 2 /(mi + m 2 ) = m 2 /(l + j 21 ) is the 
reduced mass between particles 1 and 2, 7^ = mi/mj 
is the mass ratio between particles i and j, and a = 

(TT^T)' b = a \l l21+ ri2 +{ ~' c = a721 ' Sincc 13 scalcs a11 
coordinates, they are given in terms of the smoothness of 
the interaction between particles 1 and 2 and the reduced 
mass fi\ 2 . 

The Hamiltonian (|4]) represents one particle with 
scaled mass m — 1 inside a triangular potential with 
three soft walls located at y = 0, y = ^x and y = 
— \x + L/o 3 i (Sec Fig. [T]). The smoothness of the 
walls will depend not only on the form of the potentials 
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FIG. 1. Scaled soft triangle. Dashed line represent changes of 
the wall position when the softness parameter L/azi varies. 



Vvziy), ^23(2:, y), Vz\(x, y) and the smoothness parame- 
ters (7i2, 023, 031, but also 011 a,b,c which depend only 
on the masses ratios. 

The Hamiltonian Q is quite interesting since some 
general physical situations can be observed without 
choosing a specific form for the interaction potentials: 

(i) The z dependence dissapeared due to the transla- 
tional symmetry of ((2j)- It implies the conservation of 
the total linear momentum. 

(ii) the coupling between center of mass (x) and relative 
(y) coordinates appears only in the interaction potentials 
V3\(x,y) and V2z{x,y). Thus only these potentials may 
generate a chaotic dynamics inside the billiard. When 
V31 = V23 = the Hamiltonian ^ is separable and thus 
integrable, as expected. 

(iii) Softness affects simultaneously the center of mass 
and relative coordinates since the ratios 012/023, 012/031 
multiply both variables (x,y). These quantities does not 
change the internal angles of the soft billiard. 

(iv) Changes in the masses ratio will modify the pa- 
rameters a, b and c separately, so that center of mass 
and relative coordinates will vary separately. Only these 
quantities can change the internal angles of the billiard, 
(iv) Keeping the masses ratio constant, the softness 
inside the triangle billiard is tailored by the ratios 
012/023, 012/031 and L/o-31. The hard wall limit is 
obtained by making L/a^i — > 00, i. c. the wall V31 moves 
to the right maintaining its inclination and the triangle 
size increases more and more. See dashed lines in Fig. [1] 

The Hamiltonian (0| is quite general. Any interac- 
tion potential V12, V23, V31 between particles can be used: 
Coulomb, Yukawa, Wood-Saxon, exponential etc. If the 
potential obeys the softness limits mentioned above, the 
properties (i)-(iv) should be valid for any of such inte- 
ractions. Just to exem plify this we discuss next the well 
known Toda potential [201 ■ 



IV. TODA POTENTIAL 

A. Equal masses 

We consider the three-particle Toda lattice [2(| whose 
potential is obtained by using V12 = e'-- q2 ~ qi \ V23 = 



H^ = K + eS q2 - qi) + e {q *- q2) + e^ 1 "^, 



(6) 



and corresponds to three equal masses particles moving 
on a ring. The letter e stands for equal masses. Here 
the variables (71,(72, 93 are the angles of particles related 
to the origin of the ring. In addition to the total energy 
and the linear momentum, this problem is know to have a 
third non-trivial integral of motion (For more details see 
[2ll ] ) . Using the transformation ([3]) for the particular case 
of equal masses, the following two dimensional potential 
is obtained: 



V'r 
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(7) 



The equipotential lines for this potential look like from 
a soft triangle (See Fig. (1.8) from (13). Therefore the 
analogy with the soft triangular billiard is evident. The 
additional isolated integral of motion is [22| : 

/ = 8p x (pl - 3P V ) + ( Px + V3 P y)e 2 y- 2 ^* 

-2p x e- A v + (p x -V3p y )e 2 v +2VS *, (8) 

which is not related to any obvious conservation law or 
symmetry. Potential © is a nice example of the soft 
triangular description from Eq. 

B. Unequal masses 

The unequal mass Toda problem was firstly studied 
numerically [23[ for the case of two masses interacting in 
the Hamiltonian 



K M + e~ qi + e-^-^ + e« 2 - 3, 



(9) 



where u denotes unequal masses. They showed the tran- 
sition to stochasticity when mi 7^ 777,2 and confirmed 
the intcgrability for mi = 7772. The three unequal 



masses particles was analyzed [24j for a free-end lat- 
tice Hamiltonian H^ ree = K M + e^ qi ~ q ^ + e ^- q3 \ 
They found that this problem is integrable when mi = 
e(2e - l)/(2 - e),m 2 = 2e - 1 and 7773 = 1. The pa- 
rameter e must satisfy 1/2 < e < 2. The Hamilto- 
nian H^ ree describes a scattering problem and the in- 
tcgrability cannot be verified numerically. Compared to 
the potential (O, the system H^ ree does not consider 
the interaction bettwen particles 1 and 3 because it is 
along a lattice line. See also [25[ for more about the in- 
tcgrability of Toda kind lattices. The full problem of 
the Toda-like Hamiltonian with three unequal masses 
(777,1 = TU2 = TU3 = 1) was studied before [2(| and has 
the form H T = K M + e 5 ^ 1 "^) + e A q2 ~ q ^ + e^ 1 "^. Us- 
ing the Painlcve property, they found the system to be 
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intcgrablc only for mi = ma = 6 = e = 1, which is again 
the Toda case discussed above. 

Although it is not the purpose of the present work to 
study the Toda potential, it can be written in the form 
H = K M 



Vji where 



777(91-93) 



7Tr(92-9l) 



which now includes the the softness parameters Uij be- 
tween each pair interaction. Using Eqs. ([3]) in the 
potential (fTU)) we obtain the two dimensional version 



\(ay-bx)j 



(11) 



The whole mass dependence is now inside the parameters 
from the potential Vjt- The non trivial invariant © docs 
not exist anymore for the case of unequal masses. This 
was checked by testing the condition dl u /dt = 0, where 
I u is the conservative quantity, similar to Eq. ([5} , but for 
the case of unequal masses. 

The important point to mention here is that the inte- 
grability of the system (|lip depends only on the masses 
ratios m^/mi = 7773/7711 = m^/n%2 = 1 and the softness 
ratios 0-12/0-23 = C12/C13 = 1- This generalizes the result 
from [2r| . 

In a realistic problem, like quantum dots for example, 
the material of the boundaries has impurities so that the 
interaction softness between different particles may suffer 
small changes. In such cases variations of the softness ra- 
tios may destroy the integrability. This can be very nicely 
seen in Fig.[5]where the Poincare Surface of Section (PSS) 
is shown for small variations of the softness ratios. Figure 
[2ja) is the integrable case 7712/mi = 7713/7711 = 7773/7772 
1 and CT12/CT23 = C12/C13 = 1, while Fig. [^b) we used 
C12/023 = 012/C13 = 1-2, keeping the masses ratios equal 
one. A change in the dynamics is observed and the 
chaotic motion becomes stronger. It is also possible to 
change the masses ratio keeping the softness ratios equal 
one, as can be see in Fig. [5Jc) and (d). Strong chaotic 
behavior occurs when 721 = 731 = 1.2, 732 = 1, but a not 
that large chaotic motion for 721 = 1,731 = 732 = 1.2 in 

Fig. Eld). 



V. SOFT SUITABLE POTENTIALS 
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FIG. 2. Poincare Surfaces of Section (-4.0 < y < 7.0, -10 < 
Py < 10) for E = 50 and (a) Integrable Toda case 721 = 
731 = 732 = 1.0, <Ti2/o"23 = urz/crvj, — 1.0, (b) small softness 
asymmetry 721 = 731 = 732 = 1.0, £712/0-23 = £712/0-13 = 1.2 
and the masses symmetric cases (ai2/f723 = crx%la\z = 1.0) 
(c) 721 = 731 = 1-2,732 = 1.0 and (d) 721 = 1,731 = 732 = 
1.2. 



V erf (x) = ^fErf(~ 

(X 
— 



F erf (x) = -e-^/^, 



Farc{.%) 



In the hard wall limit e — > the forces F exp (x), F er f(x) 
and F arc (x) approach the well defined S- function. For 
example, Fig. [3] shows the (a) potential V er f(x) and (b) 
force F er f{x) for distinct softness values. They show 
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To study the apropriate soft hard transition, in this 
work we assume potential functions which in limit of 
hard-walls arc not J-functions but the corresponding 
forces are. In this way the potential and forces are 
well defined for any softness of the potential. Consider, 
for example, the following three potentials and the cor- 
responding forces: 

Vexpix} — 6 ^ j F eX p{x^j — e I ^ , 



FIG. 3. (Color online) (a) Potential V er f(x) and (b) force 
F er f(x) for softness values e = 0.1 (continuous line) e = 0.05 
(dashed) and e = 0.01 (points). 

how the transition from soft to hard walls potentials and 
forces are well behaved. There are obviously other func- 
tions which in the limit e — > behave like (^-functions, but 
they are not appropriate for the purpose of the present 
work which considers billiard walls. To avoid a too long 
paper, we just discuss the potentials V exp (x) and V er f(x). 
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Results for V arc (x) arc approximately similar (not shown 
here) to V er f(x). 

Applying the potentials V exp (x) and V er f(x) in the 
Hamiltonian (jj) we obtain respectively the following sys- 
tems 



R exp = K + 2 



( — y) , — fai/- bx) . - 

d\ iJ ) I g CT 2 3 ^ y ' _L_ (j c 



= /v + Erf(-y) + £r/ 
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(ay — bx) 
+ V , 



(12) 



For -ff^ we added a constant potential Vq = 3.0 so that 
both systems have the same total bounded energy 0.0 < 
Eb < 2.0. For Eb > 2.0 the dynamics is unbounded. 



A. The soft triangles 



In this section we show that the potential energy from 
H^ p and H e ^ have the form of a soft triangle. These 
potentials are shown, respectively in Fig. @] for the case 




FIG. 5. (Color online) Potential energy from H e B xp (left) 



and H B r (right) for the case of equal softness ratio and 
721 = 5.0,731 = 1.0 and 732 = 1.0/5.0. Top row for the 
soft case Lfazx = 10 and right bottom row for Ljaz\ = 100 
approaching the hard-wall case. 



To make a distinction between the hard- wall case from 
the billiard system ([T]), here we say that when L/a^i in- 
creases, the quasi-hard wall limit is reached. In addition, 
even though we are using appropriate potentials, when 
L/a^i > 100 numerical difficults may appear because 
the potentials becomes too steep. 




FIG. 4. (Color online) Potential energy from H B xp (left) and 
H*b* (right) for the case of equal masses ratio and equal soft- 
ness ratio. Top row for the soft case L/a-n = 10 and botton 
row for L/031 = 100 approaching the hard-wall case. 

of equal softness and equal masses ratio. Top row for 
L/dsi = 10 and bottom row for L/a^i = 100. First ob- 
servation is that both potentials have a triangle billiard- 
like form (see projection on the x, y plane) and that the 
walls are soft for L/a^i = 10. As this ratio increases 
to L/a^i = 100 the walls start to look more similar to 
hard- walls. For the purpose of clarity, in the left plots 
we used the modulus in the argument of the exponential 
potentials. Otherwise it would not be possible to plot the 
coordinates in an adequate way to see the triangle form. 

For unequal masses 721 = 5.0,731 = 1.0 and 732 = 
1.0/5.0 the internal angles of the soft billiard change, as 
can be seen in Fig. [5j where the softness ratio are con- 
stant, and top row shows results for L/a^i = 10 while 
bottom row for L/a-^i = 100. Still we observe that 
L/a^i = 100 nicely approaches the hard- wall limit. 



B. The dynamics 

To analyze the dynamics related to the Hamiltoni- 
ans from (|12p we look at the PPS for the total energy 
Eb = 2.0, which is the higher energy possible before 
the particle becomes unbounded. We start showing the 
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FIG. 6. PSSs for Hg p (left) and H e B rf (right) for the case of 
equal masses and equal softness ratio. Top PSSs are the soft 
case L/1T31 = 5.0 and bottom PSSs for the quasi-hard- wall 
case L/o-31 = 20. 
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case of equal masses and equal softness ratios. Fig- 
ure [6] compares the PSSs (p y ,y,Px > 0) for the system 
H e B xp in Figs. Eta), (c), H e B rf in Figs. EIb),(d). The PSSs 
from Figs. [BJa),(b) are for the soft case Lja^x = 5.0, 
while those from Figs. [(3c), (d) for the quasi- hard- wall 
= 20. Figures Eta), (c) are related to the Toda po- 
tential discussed in section ITVAl but now for L/uzi ^ 0. 
The dynamics is regular, independent of softness param- 
eter Z//031. In this quasi-hard integrablc limit trajecto- 
ries travel along one invariant irrational torus (straight 
horizontal line) with p y > (p y < 0), but when they 
reach the boundary with a very small softness, they invert 
the momentum and return along another invariant torus. 
The softness of the walls decide how fast this momentum 
inversion occurs, i.e. for quasi-hard walls the inversion 
occurs relatively fast [Fig. Etc)]. For the hard wall limit 
from (TTJ) it occurs instantaneously but it does not change 
the torus. This is valid for these integrablc systems. 
When the dynamics from H B ^ is analyzed we see an al- 
most regular motion in the soft walls case from Fig.EJb). 
However, some hyperbolic points are present which gen- 
erate a chaotic motion when L/a^i increases. This is 
shown in Fig. Etd) for the quasi-hard case L/031 = 20, 
where the motion is chaotic with two sticky motions close 
to p y ~ 0.0, ±1.7. The sticky motion is the consequences 
of the periodic motion which occurs and can be clearly 
seen in the exponential quasi-hard limit for these points 
[Fig. Et c )]- Besides these points the motion is totally 
chaotic. This nicely shows the numerical evidence that 
the intestability condition found by the original paper 
from [20J is only valid when an exponential function is 
used. Therefore, even in the quasi- hard case L/a^i — > 00 
small changes in the potential interaction function may 
drastically change the nature of the dynamics. A chaotic 
motion was also observed (not shown) when the V arc (x) 
was used. It would be interesting if such integrability 
conditions using the exponential potentials could be ap- 
plied to find another exact solutions in quantum scatter- 
ing problems JHHl]. 

Figure [7] shows the dynamics when the masses ratio is 
changed but for equal softness ratio. The masses ratio are 
given by 721 = 5.0,731 = 1.0 and 732 = 1.0/5.0. The left 
plots are for H B xp and the right plots for H B . Figs.[7ta) 
and (b) are for the soft case L/a^i — 5.0, showing that 
the motion is almost regular for the exponential potential 
but mainly chaotic for the error potential. Figs. [?tc) and 
(d) are for the quasi-hard wall limit L/osi = 20. In both 
cases the dynamics is almost chaotic. 



VI. CONCLUSIONS 



cllcd by (5-potcntials, for which the equations of motion 
are not well defined, we propose suitable potentials to 
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FIG. 7. PSS for the case of masses ratio 721 = 5.0,731 = 1.0 
and 732 = 1.0/5.0 and equal softness ratio for H^ v (left), 
H'b (right). Above is the soft case L/a-ji = 5 and below the 
hard- wall case L/asi — 20. 



study the soft to hard interaction transition, where the 
forces, not the potentials, become 5-functions in the limit 
of hard interactions. A scaled Hamiltonian is obtained 
which nicely shows this transition and gives general clues 
about the relevant properties of interacting particles: the 
dynamics depends only on the masses ratio between par- 
ticles and softness ratio of the interaction, independent 
of the kind of interaction. Beside that any interaction 
potential with a parameter controlling the hard interact- 
ing limit can be used in this Hamiltonian which defines 
a soft triangle billiard inside which the whole dynamics 
occurs. The case of equal masses ratios is found to be in- 
tcgrable when a soft exponential interaction is assumed, 
like a Toda potential with a softness parameter. Equal 
masses ratio with other interaction potentials, like the 
error function and the arctan function (not shown) gen- 
erate a chaotic dynamics. We consider the dynamics of 
two suitable potentials, the exponential interaction with 
an additional softness parameter, a Toda like potential, 
and the error function. Results show that such potentials 
are appropriate and they do not present numerically di- 
vergencies when approaching the hard interaction limit. 



The transition from soft to hard interaction between 
three particles on a frictionlcss ring is discussed. It is 
a generalization of the two soft interacting particles dis- 
cussed in Q. Since usually hard interactions arc mod- 
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